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ABSTRACT 

Population III stars are theoretically expected to be prominent around redshifts z ~ 20, consisting of mainly 
very massive stars with M* > 100 M Q , but there is no direct observational evidence for these objects. They 
may produce collapsar gamma-ray bursts (GRBs), with jets driven by magnetohydrodynamic processes, whose 
£NJ , total isotropic-equivalent energy could be as high as > 10 57 erg over a cosmological -rest-frame duration of 

^ ■ td > 10 4 s, depending on the progenitor mass. Here we calculate the afterglow spectra of such Pop. Ill GRBs 

based on the standard external shock model, and show that they will be detectable with the Swift BAT/XRT and 
Fermi LAT instruments. We find that in some cases a spectral break due to electron-positron pair creation will 
be observable in the LAT energy range, which can put constraints on the ambient density of the pre-collapse 
Pop. Ill star. Thus, high redshift GRB afterglow observations could be unique and powerful probes of the 
properties of Pop. Ill stars and their environments. We examine the trigger threshold of the BAT instrument in 
detail, focusing on the image trigger system, and show that the prompt emission of Pop. Ill GRBs could also 
be detected by BAT. Finally we briefly show that the late-time radio afterglows of Pop. Ill GRBs for typical 
parameters, despite the large distances, can be very bright: ~ 140 mJy at 1 GHz, which may lead to a constraint 
on the Pop. Ill GRB rate from the current radio survey data, and ~ 2.4 mJy at 70 MHz, which implies that 
Pop. Ill GRB radio afterglows could be interesting background source candidates for 21 cm absorption line 
O . detections. 

Subject headings: black hole physics — dark ages, reionization, first stars — gamma rays burst: general — 
stars: Population III — X-rays: bursts — radio continuum: general — surveys 

a 

1. INTRODUCTION 

\ Recent studies on cosmology and primordial star formation predict that the first generation of stars (population III stars) may 
be most prominent around z ~ 20, consisting of metal-poor, mainly very massive stars (VMSs) with M* > 100 Mq (e.g., Abel et 
Q\ * al. 2002; Omukai & Palla 2003; Yoshida et al. 2006; Ciardi & Ferrara 2005). These first stars are thought to play a significant 
\Q , role in setting off cosmic reionization, in the initial enrichment of the intergalactic medium (IGM) with heavy elements, and in 
C\J ■ seeding the intermediate and supermassive black holes (BHs) encountered in galaxies. The details of how these processes unfold 
t-H ' remain elusive, since observational data for redshifts z > 6 are very limited. 

Observations of gamma-ray bursts (GRBs), however, may provide unique probes of the physical conditions of the universe at 
• such redshifts. The GRB prompt emission and the afterglows were expected to be observable at least out to z > 10, with their 
' redshifts being determined through the detection of a Lya drop-off in the infrared (IR), or through redshifted atomic lines back- 
i— ( , lighted by the afterglows (e.g., Lamb & Reichart 2000; Ciardi & Loeb 2000; Gou et al. 2004). This can serve as a tracer of the 
!• ■ history of the cosmic star formation rate (e.g., Totani 1997; Porciani & Madau 2001; Bromm & Loeb 2006; Kistler et al. 2009), 
.£h providing invaluable information about the physical conditions in the IGM of the very high redshift universe (e.g., Barkana & 
, Loeb 2004; Ioka & Meszaros 2005; Inoue et al. 2006). Currently the most distant object that has been spectroscopically confirmed 
is GRB 090423 at z ~ 8.2 (Tanvir et al. 2009; Salvaterra et al. 2009), and the detailed spectroscopic observation of GRB 050904 
d ' at z — 6.3 has put a unique upper bound on the neutral hydrogen fraction in the IGM at that redshift (Totani et al. 2006; Kawai 
et al. 2006), indicating that GRB observations are very promising for exploring the high-redshift universe (see also Greiner et al. 
2009, for GRB 080913 with z ~ 6.7). 

In those previous papers, the GRBs arising from Population III VMSs were considered to have similar properties as the 
GRBs arising from Population I/II stars, e.g., they were usually assumed to have similar luminosity functions, even if perhaps 
extending to somewhat higher masses, and most importantly, their radiation properties, durations and spectra were modeled as 
being essentially similar to their lower redshift counterparts. However, this simplifying assumption may not be valid, as pointed 
out by Fryer et al. (2001) and Komissarov & Barkov (2010) (hereafter KB 10). One difference is that the accretion disks around 
the much larger black holes resulting from core collapse of VMS progenitor would be too cool to lead to neutrino-cooled thin 
disks and conversion of neutrinos into electron-positron pairs (Eichler et al. 1989; Woosley 1993). Thus, the Pop. Ill GRBs are 
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much likelier to be driven by MHD processes, converting the rotational energy of the central BH into a Poynting-flux-dominated 
jet (Blandford & Znajek 1977), rather than the usually assumed thermal-energy-dominated jets. The total energy of a Pop. Ill 
GRB is then proportional to the total disk (torus) mass, which in turn can be assumed to be proportional to the progenitor stellar 
mass, which can be much higher than that of a Pop. I/II GRB. In addition, the fall-back time and/or the disk accretion time, i.e., 
the active duration time of a Pop. Ill GRB jet can be much longer than that of a Pop. I/II GRB jet, due to the larger progenitor 
star. 

Building on this premise, Meszaros & Rees (2010) (hereafter MR10) proposed a possible model of the prompt emission and 
afterglow of such Pop. Ill GRBs, and made rough predictions for their observational properties in the Swift and Fermi satellite 
bands. In this paper, we calculate in significantly more detail the very early afterglow properties of Pop. Ill GRBs, and show 
that the combination of Swift and Fermi observations, complemented by deep IR observations of the afterglow immediately 
following the prompt emission, can constrain the total isotropic-equivalent energies (Ei S0 ) of the Pop. Ill GRBs, as well as the 
particle densities n of their circumburst medium. The detection of a burst with a very high total isotropic-equivalent energy 
Eiso ^ 10 57 erg and a very long (cosmological rest frame) duration td > 10 4 s would be strong evidence for a VMS progenitor. To 
constrain the total energy and the duration, observations of the prompt emission, whose interpretation is more dependent on the 
model details, should be complemented with observations of the afterglow, which is much less model-dependent. 

The properties of the circumburst medium, i.e., the environments of the first stars prior to their collapse, have so far only been 
inferred from model numerical simulations, which differ significantly among each other. For example, the typical galactic gas 
environment could evolve as n oc (1+z) 4 (Ciardi & Loeb 2000), or it might be approximately independent of redshift, n ~ 0.1 cm" 3 , 
as a result of stellar radiation feedback (Whalen et al. 2004; Alvarez et al. 2006). Observations and modeling of high-redshift 
GRB afterglows could distinguish between such numerical models. The small number of analyses of what are currently the 
most distant GRBs imply that the circumburst densities of these high-redshift GRBs could be very different from each other, e.g. 
n ~ 10 2 - 10 3 cm" 3 for GRB 050904 with z ~ 6.3 (Gou et al. 2007), and n ~ 1 cm" 3 for GRB 090423 with z ~ 8.2 (Chandra 
et al. 2010). A high total energy and a high circumburst medium density could lead, in principle, to such a high compactness 
parameter of the shocked afterglow emission region that a spectral break due to e + e~ pair production (i.e., 77 self-absorption 
break) may be observable in the Fermi LAT energy range. This is a unique and interesting point, which we explore here, since 
the 77 self-absorption is usually not significant for the afterglows of Pop. I/II GRBs (Zhang & Meszaros 2001). Using this, we 
could constrain the circumburst density n from the observation of a 77 self-absorption break, which is a new method to constrain 
the environment of the pre-explosion GRB host galaxy. 

The external shock model of the GRB afterglows seems to be robust, since it can explain many of the late-time multi-band 
afterglows detected so far, and a simple extension of this model (e.g., continuous energy injection into the external shock) may 
explain many of the early-time (observer's time f b s < 1 hr) X-ray and optical afterglows detected by Swift (Liang et al. 2007). 
More importantly, some of the very early high-energy afterglows (immediately following the prompt emission) recently observed 
by Fermi LAT are shown to be explained by this model (e.g., Kumar & Barniol Duran 2009; De Pasquale et al. 2010; Corsi et al. 
2010). 

The basic parameters of the Poynting-dominated Pop. Ill GRB model are defined in Section 2. The very early afterglow 
spectrum is calculated in Section 3 (based on the standard external shock model described in Appendix), where we discuss 
how to constrain Ei S0 ,td and n from the observations. In Section 4 we deduce the effective trigger threshold of the Swift BAT 
instrument, and show that the prompt emission of Pop. Ill GRBs could trigger BAT. In Section 5 we compute the late-time radio 
afterglow flux for a typical set of parameters, and evaluate the current radio survey data constraints on the Pop. Ill GRB rate, as 
well as the prospects for 21 cm absorption line detection in the Pop. Ill GRB radio afterglow spectra. In Section 6 we present a 
summary of our findings. 

2. POYNTING-DOMINATED POP. Ill GRB MODEL 

We consider VMSs rotating very fast, close to the break-up speed, as a representative case of Pop. Ill GRB progenitor stars. 
Those in the 140 M Q < < 260 M range are expected to explode as pair instability supernovae without leaving any compact 
remnant behind, while those in the M* > 260 M range are expected to undergo a core collapse leading directly to a central 
BH, whose mass would itself be hundreds of solar masses (Fryer et al. 2001; Heger et al. 2003; Ohkubo et al. 2006). Accretion 
onto such BHs could lead to collapsar GRBs (Woosley 1993; MacFadyen & Woosley 1999). Prior to the collapse, the fast 
rotating VMSs may be chemically homogeneous and compact, without entering the red giant phase, so that the stellar radius is 
/?* ~ 10 12 cm for M* ~ 10 3 M Q (KB10; Yoon et al. 2006; Woosley & Heger 2006). 

For such large BH masses M/, > 100M Q , the density and temperature of the accretion disk are too low for neutrino cooling to be 
important, and the low neutrino release from the accretion disk is insufficient to power a strong jet (Fryer et al. 2001). The rate of 
energy deposition through this mechanism may be estimated by using the formula recently deduced by Zalamea & Beloborodov 
(2010) 

L„ p ~ 5 x 1 46 M% A M~^l erg s" 1 , ( 1 ) 

where M = Q.\M-\Mq s -1 is the accretion rate and Mh,2.5 = Mi,/(10 25 Mq). This is clearly insufficient for detection from such 
high redshifts. However, strong magnetic field build-up in the accretion torus or disk could lead to much stronger jets, dominated 
by Poynting flux. Such jets will be highly relativistic, driven by the magnetic extraction of the rotational energy of the central 
BH through the Blandford-Znajek (BZ) mechanism (Blandford & Znajek 1977). The luminosity extracted from a Kerr BH with 
dimensionless spin parameter a/, threaded by a magnetic field of strength £/, is (Thorne et al. 1986) 

L BZ ^^-B 2 h R 2 h c, (2) 



Population III GRB Afterglows 



3 



Prompt Emission 






External Shock 




\ Emi s sion 











t d (1+z) -lday log(t obs ) 



FIG. 1 . — Schematic figure of the bolometric luminosity evolutions of the external shock emission (thick line) and the prompt emission (thin line) as functions 
of time from the onset of the prompt emission. The prompt emission may involve some variability, which is shown by the dotted line. In this figure we assume 
that the jet luminosity Lbz is constant for t ^ s < f,; obs = +z) an d subsequently declines rapidly. The external shock is in the thick-shell regime. Its luminosity 
is also constant for t y, s < tj „\, s and subsequently enters the self-similar phase L <x t~^ s . The argument is similar for the case of the jet luminosity Lbz r q with 
q < 1, where luminosities of the prompt emission and external shock emission both evolve as L oc before tj,^ and in the same way as this figure after t^obs- 



where Rh ~ GMh /c 2 is the event horizon radius of the BH. The dynamics of the radiatively inefficient accretion disk may be 
described (KB10) through advection-dominated (ADAF) model (Narayan & Yi 1994). For a VMS rotating at, say, half the 
break-up speed, the disk outer radius will be R d ~ R*/4, and for a disk viscosity parameter a = lCT'cu-i, the accretion time is 

^3^(^) 1/ ^ 1 - 4Xl04a= ^ M ^ S ' (3) 

where we have defined R* : n = R*/i0 12 cm. This gives an estimate for both the disk lifetime and the duration of the jet, in 
the source frame. Given the jet propagation speed inside the star, <~ 0.2c, deduced from magnetohydrodynamic simulations 
(Barkov & Komissarov 2008), the intrinsic jet duration t d ~ 10 s is sufficient to break through the star. The poloidal magnetic 
field strength in the disk should scale with the disk gas pressure, P, so that B\ = 8irP/f3, where j3 = 10 /3\ is the magnetization 
parameter (e.g., Reynolds et al. 2006). Then we have 

Combining these equations result in a jet luminosity 



L BZ * ^-^Mc 2 ~ 2.2 x 10 51 -^M d , 2 . 5 t d \ erg s" 1 , (5) 

where for the second equalities in Eqs. (4) and (5) we have assumed a constant accretion rate M ~ M d / td, and Md = IO^M^sMq 
is the total disk mass. 

Let us assume that the factor af,/(aP) is roughly constant, so that Lbz oc M. The total extracted energy during the accretion 

time td is then L B z — (\/l4/96)(a|/ a/3)MdC 2 . Assuming that the jet has an opening angle of 9j = 0.1 we can then write the 
total isotropic-equivalent energy of the jet as 

L IS0 ~ 4.4 x 10 57 <±Z±g!L M fa erg, (6) 
a-ipi 

where e 7 is the radiation efficiency of the prompt emission, and 1 - e 7 is of order of unity. Equation (6) is also applicable to cases 
where M is not constant. As far as L B z oc t~ q (i.e., M oc t~ q ) with q < 1, the jet can break out from the star and subsequently keep 
injecting energy into the external medium for the duration of the order of td- The forward shock produced in the external medium 
enters a self-similar expansion phase with total shock energy ~ Li S0 soon after t = td (Blandford & McKee 1976). Interestingly, 
the value of L; so for a disk mass Md ~ 3 M is consistent with the observed largest value of the isotropic-equivalent 7-ray energy 
£ 7iiso ~ 10 55 erg for GRB 080916C at the redshift z ~ 4.4 (Abdo et al. 2009). (On the other hand, the isotropic luminosity 
Lj so ~ 4.4 x 10 53 {a 2 h / af3)Md2.5Q\0~ 2 _ l erg s" 1 is comparable to the observed largest value.) Thus, if we were to observe a burst 

at redshift z > 10 with E iso > 10 57 erg, and with a self-similar phase starting at t d , obs = t d (l+z)> 1 day, this would very likely be 
a burst from a Pop. Ill VMS with > 300 M Q . 
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FIG. 2. — Example of the observer-frame spectrum of a Pop. Ill GRB at the time tj b s when the jet activity ends, for the case of negligible intra-source pair 
production. The parameters are £57. 6 = tj 4 = no = £b-2 = £e.-l = f(.P) = 1, and the source redshift 1 +z = 20. The dot-dashed line shows the external shock 
(afterglow) spectrum, which consists of the synchrotron and SSC components (solid lines). The synchrotron component peaks at e m ~ 5.4 eV, and the SSC 
component peaks at e^ c ~ 1.8 X 10 2 MeV. The 77 self-absorption break energy is £ 77 ~ 21 GeV, which is larger than the SSC peak eff, so that most of the SSC 
emission escapes without being absorbed within the emitting region. The 77 absorption due to the EBL is expected to become significant at e > £ebl — 7 GeV 
(Inoue et al. 2010), as shown by the dashed line. The dotted line represents the prompt emission's dominant photospheric black-body component with a possible 
power-law extension, assuming l+a = 10 and L 53 6 = r; 8 = T; = 1 (see Section 4.2). 



3. VERY EARLY AFTERGLOW SPECTRUM 

Rough predictions for the observational properties of the afterglows of Pop. Ill GRBs were made in MR10. The external shock 
driven by the jet in the circumburst medium powers the afterglow, which can be studied independently of the prompt emission 
(Meszaros & Rees 1997a; Sari et al. 1998; Sari & Esin 2001). This is true whether the jet is baryonic or Poynting-dominated, 
the jet acting simply as a piston. 7 The external shock amplifies the magnetic field in the shocked region via plasma and/or 
magnetohydrodynamic instabilities, and accelerates the electrons in the shocked region to a power-law energy distribution. The 
accelerated electrons produce synchrotron and synchrotron-self-Compton (SSC) radiation as an afterglow. Here we go beyond 
the previous schematic outlines, and calculate the spectrum of this emission in detail, including the Klein-Nishina as well as e + e~ 
pair formation effects. 

The bolometric luminosity of the external shock emission (with prompt emission light curve) is illustrated in Figure 1 . We 
focus on the external shock emission at the observer's time f D b s ~ f^obs, near the beginning of the self-similar expansion phase of 
the shock, when the emission is bright and may not be hidden by the prompt emission. 

Calculations of the external shock emission spectrum involve the parameters E lso and fj, as well as the external medium 
number density n, the fractions e# and e e of the thermal energy in the shocked region that are carried by the magnetic field and 
the electrons, respectively, and the index p of the energy spectrum of the accelerated electrons. We calculate the external shock 
emission spectrum of a Pop. Ill GRB at f^bs based on the standard model described in Appendix. 

As introduced in Section 1, the circumburst medium density in the very high-redshift universe is likely to be n > 0.1 cm" 3 . The 
microphysical parameters maybe independent of E vio ,td 1 orn as long as the shock velocity is highly relativistic, so that e#,e c , and 
p are thought to be similar to those for the bursts observed so far. Those have been constrained by fitting the late-time afterglows 
through models (which are similar to our model shown in Appendix). The parameters related to the electrons are constrained 
relatively tightly as e e ~ 0.1 and p <~ 2.3, while those for the magnetic field are not so tightly constrained, although typically for 
many afterglows 10~ 3 < eg < 10 _1 (e.g., Panaitescu & Kumar 2002; Wijers & Galama 1999). 

The external shock emission at Q b s will have two intrinsically different cases, depending on the significance of the electron- 
positron pair creation within the emitting region. These cases are characterized by a negligible pair production regime and a 
significant pair production regime, which we show examples of spectra separately below. 

3.1. Case of Negligible Pair Production 
An example of the negligible pair production case is obtained for the parameters 

£57.6 = t d A = «o = £fl -2 = e e -1 = f(p) = 1 , (7) 

7 The same is not true for a reverse shock, whose existence and properties are more dependent on the nature of the ejecta jet (Mimica et al. 2009; Mizuno et al. 
2009; Lyutikov 2010). However, a reverse shock emission is most prominent in the low frequencies, e.g., the IR bands, while we are interested in the X-ray and 
7-ray bands at tj, so this is not considered here. 
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FIG. 3. — Example of a Pop. Ill GRB spectrum at time td with significant pair production in the external shock. The parameters are £57.6 = *rf.4 = €B-2=f(p) = 1. 
n = 10 2 , and e e _i = 2, and the source redshift 1 +z = 20. The dot-dashed line shows the external shock emission spectrum, which consists of the four components 
(solid lines). The synchrotron component of the original electrons peaks at e m ~ 21 eV. The peak of the SSC emission of the original electrons (not shown in 
this figure) ejjf ~ 9.2 X 10 2 MeV is above the 77 self-absorption break energy e~ n ~ 6.1 X 10 2 MeV, so that most of the SSC emission is absorbed within 
the emission region. The pairs emit synchrotron and SSC radiations peaking at e±. p ~ 0.11 eV and e^r p ~ 23 keV, respectively. The IC-scattered original 
synchrotron emission by the pairs and the IC-scattered pair synchrotron emission by the original electrons have spectra with similar characteristic energies 



(including the peak energies at e„, t 



Z ±,P - 



■ 4.5 MeV) and different flux normalizations, which have been superposed in this figure. The 77 absorption due to 



the EBL is expected to be significant above ~ 7 GeV (Inoue et al. 2010), as shown by the dashed line. The dotted line represents the prompt emission's dominant 
photospheric black-body component with a possible power-law extension, assuming 1 + cr = 10 and L53 6 = r/ g = T; = 1 (see Section 4.2). 



where the notation Q = 10*2* i n c g s units has been adopted (£57.6 = £i SO /T0 576 erg). The overall observer-frame spectrum for 
this case is shown in Figure 2 (see Appendix A.l). The synchrotron emission spectrum peaks at e m ~ 5.4 eV with the flux 
e m F £m ~ 4.4 x 10" 10 erg cm" 2 s" 1 , having spectral breaks at synchrotron self-absorption energy e a ~ 9.7 x 10" 4 eV and at energy 
corresponding to the maximum electron energy Em — 1 -8 x 10 2 MeV. The SSC emission spectrum peaks at e^ c ~ 1 .8 x 10 2 MeV 
with the flux ef^F^ ~ 1.2 x 10" 9 erg cm" 2 s" 1 and has a spectral break at 77 self-absorption energy e 77 ~ 21 GeV. Since 
£ 77 > 4°, most of the SSC emission is observed without being converted into e + e~ within the emitting region. 

Even photons escaping without attenuation within the emitting region can be absorbed by interacting with the extragalac- 
tic background light (EBL) (77 — > e + e~). Inoue et al. (2010) use a semi-analytic model of the evolving EBL and expect that 
high-energy photon absorption by the EBL for an arbitrary source at z — 20 is significant at e > eebl — 7 GeV. For the param- 
eters adopted here, e 77 is larger than Eebl, which precludes obtaining intrinsic information about the emitting region from the 
observation of the 77 break. 

The temporal evolution of the characteristic quantities during the self-similar phase, i.e., at f t, s > ^,obs is obtained replacing 
td by the variable f b s in the model equations in Appendix, and taking all other parameters as constant. On the other hand, for 
fobs < fj.obs, if the jet luminosity evolves as Lbz oc t~ q with q < 1, one can obtain the temporal evolution of the characteristic 
quantities by taking £i S0 oc t^, replacing td by f bs> an d taking all the other parameters as constant. In the high-energy range, 



e?, c < e < e 



77 



as an example, we obtain F 



sc 



oc t. 



(3/4)(2-p)-(3/8) 9 (p+2/3) 
obs 



for f obs < t dfi b s , and oc t ( J^~ 9p)/8 for f obs > t d , bs, which 



implies a steepening break at fj j0 b s for p ~ 2 and q < 1. Thus one can identify the jet duration f^ obs as the observed break time. 
(One can also estimate f^ obs by the duration of the prompt emission.) 

3.2. Case of Significant Pair Production 
An example of the case of significant pair production is obtained with the parameter set 

£ , 57.6 = ^,4 = eB-2=/0) = l, «o = 10 2 , e e _i=2. (8) 

The overall observer-frame spectrum for this case is shown in Figure 3 (see Appendix A. 2). The synchrotron emission spectrum 
of the original electrons peaks at e m ~ 21 eV with the flux e m F £m ~ 5.8 x 10" 10 erg cm" 2 s" 1 , having the maximum energy 
Em — 1 -0 x 10 2 MeV. The 77 self-absorption energy is e 77 ~ 6. 1 x 10 2 MeV, and thus most of the SSC emission with the spectral 
peak 4 C ~ 9.2 x 10 2 MeV with the peak flux 4 C F £ SC - 2.7 x 10" 9 erg cm 2 s 1 is absorbed within the emitting region. The 
created pairs emit synchrotron emission peaking at e±_ p ~ 0.1 1 eV and SSC emission peaking at e^ p ~ 23 keV. In addition to 
these, the pairs Inverse Compton (IC)-scatter the original electron synchrotron emission and the original electrons IC-scatter the 



pair synchrotron emission, which have similar spectra peaking at e± p ■ 
been superposed in Figure 3. 



-IC , 



: 4.5 MeV with different flux values, and they have 
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FIG. 4. — The spectra for a lower redshift 1 +z = 10, showing the total external shock emission in the case of negligible pairs (solid curve; c.f. Figure 2) and in 
the case of significant pairs (dot-dashed curve; c.f. Figure 3), with the prompt photospheric emission (dotted curve). The EBL 77 cutoff is shown by the dashed 
lines, with a cutoff energy eebl assumed to be similar to that in the 1 + z = 20 figures. 



In this case, we can measure e 77 , since this is well below the expected EBL cut-off, so one would be able to draw inferences 
about the source parameters from the 77 self-absorption break (see next section for details). 
The temporal evolution of the flux in the LAT energy range will be F e ic <e<e _ ri oc t ^J %){p 2) {l l x6 ^ lp+T> f or t ohs < i d obs anc [ ^ 

t das~ 21P ^ 16 f° rf <*s > f^obs, which implies a steepening break at fj, b s . The lightcurve well after td,obs, however, may be complicated. 
The SSC component could become dominant at later times, since the SSC energy evolves as e^ c ex tj^ 4 and the 77 self-absorption 
energy evolves as e 77 oc tf b { 4) ~ (2/p) for f obs > t d:0bs . 

3.3. Constraints on Ei S0 and n from Observations 

We have shown two typical cases of the external shock emission of Pop. Ill GRBs, one being the case of negligible pair pro- 
duction, and the other being the case of significant pair production for which the cascade process stops when the first generation 
pairs are created. There may be cases where the cascade process can create second (or higher) generation pairs, e.g., for larger e e 
and/or larger external n (in the above example we used e e = 0.2 and n = 10 2 cm" 3 ). In any case, the important point is that one will 
be able to detect a spectral break at energy e 77 due to pair creation within the emission region in the Fermi LAT energy range, 
50 MeV - 30 GeV. This is a unique feature of GRB afterglows with very large £i S0 (as expected for Pop. Ill GRB) and modest 
to large external density n. Equation (A18) or (A22) indicate that larger £i S0 and n lead to smaller e 77 , increasing its diagnostic 
value. This is in contrast to the usual case of Pop. I/II GRBs, where the 77 self-absorption energy is not relevant for observations 
(Zhang & Meszaros 2001). 

We can estimate the detection thresholds in the high energy ranges from the joint observation of GRB 090510 by Swift and 
Fermi (De Pasquale et al. 2010). This indicates that the thresholds of the 1-day averaged sF e flux are <~ 6 x 10~ 15 erg cm" 2 s" 1 in 
the XRT energy range 0.3 - 10 keV, ~3x 10~ 10 erg cm" 2 s" 1 in the BAT energy range 15 - 150 keV, and <~ 3 x 10" 1 1 erg cm" 2 s" 1 
in the LAT energy range 50 MeV - 30 GeV. Compared to the results shown in Figures 2 and 3, it appears that the thresholds of the 
XRT and LAT are thus sufficiently low, and the BAT is marginally low only for the case of Figure 3, to observe the high-energy 
spectrum of the external shock emission of Pop. Ill GRB. Furthermore, for both cases of Figures 2 and 3, we find that the 
very high energy emission at e > 10 GeV could be detected with next generation facility such as Cherenkov Telescope Array 
(CTA) 8 , which will have a threshold of 1-day averaged eF £ flux ~ 10~ 12 erg cm" 2 s" 1 for e > 10 GeV, although this significantly 
depends on the EBL attenuation for individual burst. In Figure 4 we show, for reference, the similar results for a lower redshift of 
1 +z = 10, using the same values of the other parameters as for the previous two figures. The above statement is also applicable to 
this case. The EBL attenuation effect for 1 +z = 10 is expected to be similar to that for 1 + z = 20 since the EBL intensity declines 
at z > 10 (Inoue et al. 2010). 

One of the main questions that will be asked, if and when the redshift of a burst is determined to be z > 10 (e.g. by observation 
of the Lyct cutoff at IR frequencies), is whether this burst is produced by a Pop. Ill VMS or not. An effective way to pinpoint 
a Pop. Ill progenitor is examine whether the afterglow spectrum from its surrounding medium is devoid of metals through high 
resolution IR and X-ray spectroscopy by ground based facilities and/or future space experiments. Here we detailedly discuss 
another way by estimating the duration and total energy of the jet through the X-ray and 7-ray observations of the afterglow 

8 http://www.cta-observatory.org. 
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and/or prompt emission. The jet duration timescale td can be estimated from the steepening of the afterglow light curve and/or 
the end of the prompt emission (see Sections 3.1 and 3.2, and Figure 1). A lower bound on the total isotropic-equivalent energy 
iiiso can be estimated from the observed flux level in a specific energy range. Thus, if one obtains an £i S0 > 10 57 erg and td > 10 4 s, 
this burst is almost certainly bound to be a Pop. Ill GRB. 

One could in principle attempt to constrain the electron spectral index p as well, using the observed photon spectral index, and 
the four physical parameters E{ S0 , n, e B , and e e , by measuring the observables e m , F £m , e^ c and with the Swift XRT, the Fermi 
LAT, and a possible IR detection. However, we have shown that there are some cases in which the spectral peak in the LAT energy 
range corresponds to £± p . In such cases one cannot easily distinguish between the regime of negligible pair production and that 
of significant pair production, so that one would not be able to uniquely constrain all four physical parameters. Furthermore, the 
prompt emission may last until f D b s ~ f^obs, which could hide the peak of the external shock synchrotron or IC/SSC emission, 
reducing the number of observables of the external shock emission. 

We consider now the case in which the external shock emission is observed in the LAT energy range, without being hidden by 
the prompt emission component, and show how this can constrain £j so as well as n. In this case we have two observables, the flux 
at some energy in the LAT energy range, SlF El , and the 77 self-absorption break energy e 77 , which is identifiable if it is well 
below £ EB l- 

The LAT flux can put a lower limit on E iso , modulo the uncertainty on e e . The LAT flux should be lower than e^F^ or 
£ ±,p^e± p f° r tne case °f e « ^ e s or e m F Sm for the case of e e <C £b- In any case, we have that SlF Sl < e e Ei so t d ~ 1 (p-2)/[4irdl(p- 1)]. 
We may approximate p to be 2X L , where A/, is the measured spectral index in the LAT energy range. This leads to 

£ lso > 2 x 10 57 e~U ( 10 ^,^-2 s -i ) ^/l(A l )4 20 erg, (9) 

where /l(Al) = (2/7)(Al— 0.5)/(Al- 1), and the luminosity distance cIl is normalized by the value for 1 +z = 20, 6.7 x 10 29 cm. 
This bound can be compared with the total isotropic-equivalent energy of the prompt emission, E^^. 

From an observation of the 77 self-absorption break in the LAT band, we can then constrain the density of the medium around 
the Pop. Ill star before its collapse. The EBL cutoff energy £ EB l can be estimated by some EBL models (e.g., Inoue et al. 
2010) when we have the source redshift z. If we detect a spectral break which is well below the values of Eebl expected for 
practically all EBL models, that is likely to be e 77 . By using this, we can constrain the bulk Lorentz factor of the emitting 
region Td (see Appendix and Zhang & Meszaros 2001; Lithwick & Sari 2001). For the photons at e 77 in the LAT energy range, 
the main target photons have energies at e at ~ T 2 d m 2 c 4 /[(1 +z) 2 e 77 ] ~ 7 (r rf /10 2 ) 2 (e 77 /l GeV) _1 [(l +z)/20]" 2 keV. Thus we 
can estimate the target photon number density by using the Swift XRT data. The equation for the optical depth r 77 (e 77 ) ~ 
(<T T /W)(dl/ r 2 d )F Vx (e at /s x y x td/[h(\+z)} = 1, where 

£x, F Vx and A are the observed X-ray energy, flux and spectral index, with 
the equation for the emission radius ~ cT d td, leads to 

r^60(2o,W (^) ± ( T ^)»( T i^)»# < ,A [(1+z , /2 „ 1 «a. 

By using E iso ~ 4irr d Y 2 d nm p c 2 , we have an estimate of E i>i0 /n. Combining it with the above lower limit on -Eiso-. we can put a 
lower limit on n, 

. >40 (20,™ ^ ( 1 ^)*( T gv)*( I ^)* {Wfg^ , 

x tj* f L (\ L )d^[{\+z)/20]^ cm-s. (11) 

Even if the prompt emission hides the external shock X-ray emission, taking the prompt emission X-ray flux as F vx may provide 
us with a good estimate for a lower limit on n. 

In principle, one could have a situation where e e <C £b, although this appears to be rare. In this case the SSC component 
would be dim, and any high-energy cutoff (or break) is due to a synchrotron maximum energy e M . Since e 77 oc f"^ 1 for p <~ 2 

and Em oc t^, we can distinguish the cutoff origins. If the cutoff is Em, we can compute the bulk Lorentz factor by using 
Equation (A 14) and a lower limit on n. 

Above we have considered the cases in which the afterglow emission can be well observed in the XRT and LAT energy ranges. 
It would be useful to examine for what ranges of parameters the emission cannot be well observed. At e > e m , we have a rough 

but simple estimate of the flux as eF £ <~ e e E l!io t d l (p — 2) / '[4nd 2 (p - 1)] oc e e Li so dl 2 f(p), where Lj so = E[ so t d l oc 
and the flux is found to be weakly dependent on eg or n. Although we have still several free parameters, we may examine the 
rough parameter dependence of the flux level. As an example of smaller progenitor mass M» ^100 M Q (which corresponds to 
Md,2.5 — Mh,2.5 — 0.1), if we assume R*.\2 — 0.5 and Oj-i ~ e e _i ~ f(p) ~ di,20 — L we have eF £ ~ 1 x 10~ 10 erg cm -2 s" 1 at 
frf.obs ~ 3 day, which is still well above the XRT threshold, but below the BAT threshold, for the 1-day integration of the flux 
(shown above). This flux is marginally above the LAT threshold, but the spectral break at e 77 may be difficult to be clearly 
identified. However, note that the flux could be much higher than this value, depending on the poorly constrained parameters 
6j, and e e . The XRT instrument seems to be very powerful to observe the emission for even smaller mass progenitors, but 
anyway in order to identify the direction to the GRBs on the sky, the emission flux of the afterglow or the prompt emission has to 
be sufficiently high to trigger the BAT instrument. This issue is discussed below. 
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FIG. 5. — The 1-s peak photon flux in the 15-150 keV band versus the fluence in the 15-150 keV band for the Swift GRBs in the BAT2 catalog (Sakamoto et 
al. 2010) overlaid with 64s image triggered GRBs (red). 



4. DETECTABILITY AND RATE OF POP. Ill GRBS 

4.1. Swift BAT detection threshold 

Instruments with large field of view need to be triggered by a GRB in order to observe them from early times. Here we 
investigate the detection threshold of the Swift BAT in detail to address whether Pop. Ill GRBs can trigger the BAT, by focusing 
on the "image trigger" mode. 

The regular "rate trigger" mode looks for a rate increase in the light curves. The rate trigger is thus sensitive to a burst which 
is variable on a relatively short time scale. On the other hand, the image trigger searches for a burst by creating sky images every 
64 s in the 15-50keV band. The image trigger is purely based on whether a new source is found in the sky image or not in the 
given interval, without looking for a rate increase in the light curves. Generally it is not simple to find the detection threshold of 
the BAT instrument, because 494 different trigger criteria (e.g. different energy bands, time scales and combinations of detectors) 
have been running on-board. We are able to estimate a reasonable detection threshold, however, by investigating only the BAT 
image triggered GRBs, since the image trigger with the 64 s integration is based on a criterion with a fixed time-scale and energy 
band. 

There are some other reasons for focusing on the image triggered bursts to obtain a reasonable detection threshold of BAT. 
Figure 5 shows a diagram of the peak photon fluxes and the fluences in the 15- 150 keV band of Swift GRBs in the BAT2 
catalog which includes 64 s image triggered GRBs (Sakamoto et al. 2010). The BAT GRBs found by the 64 s image trigger have 
systematically weaker peak photon flux (but similar fluence) compared to those GRBs found by the rate triggers. Furthermore, 
the prompt emission and external shock emission of Pop. Ill GRBs are likely to have a lower peak flux and a longer duration, and 
to be less variable due to intrinsic property and high-redshift time stretching. The image trigger mode may be more sensitive to 
such emission than the rate trigger mode. Indeed, the high-redshift GRB 050904 (with z — 6.3, Kawai et al. 2006) and the low- 
luminosity GRB 060218 (Campana et al. 2006; Toma et al. 2007) were image triggered bursts. Thus the estimate of the threshold 
of the image trigger mode may be relevant especially for Pop. Ill GRBs. (Note, however, that GRB 090423 (with z — 8.2, Tanvir 
et al. 2009; Salvaterra et al. 2009) and GRB 080913 (with z ~ 6.7, Greiner et al. 2009) are the rate triggered bursts.) 

There are 72 GRBs (out of 467 GRBs) detected by the image trigger in the 2nd BAT GRB catalog. Out of 72 imaged triggered 
GRBs, there are 50 GRBs found by the image trigger with an integration time of 64 s, 9 which have been shown in the red points 
in Figure 5. We focus on this sample of 50 GRBs to estimate the BAT detection threshold. To understand the BAT detection 
threshold for the image trigger GRBs, we create the spectrum using the 64 s image trigger interval, and extract the photon 
fluence of the interval in the 15-50 keV band (same energy band of the image trigger) and the best fit photon index based on a 
simple power-law model. The image trigger interval has a duration of 64 s starting from the BAT trigger time. Figure 6 shows the 
distribution of the photon index and the photon fluence in the 15-50 keV band for a 64 s image interval of 50 image trigger GRBs. 
The imaged fluence in the 15-50keV band does not have a strong dependence on the observed spectrum. Therefore, we conclude 
that the photon fluence threshold of the BAT 64 s image trigger is ~ 1 ph cm" 2 in the 15-50 ke V band, corresponding to a photon 
fluence threshold in the 15-150 keV band, of ~ 2 ph cm . This corresponds to the averaged photon flux ~ 0.03 ph cm s" 1 in 
the 15-150 keV band. Note that this is the minimum averaged flux in the first 64 s interval of the image triggered GRBs, which 



Although the image trigger is basically producing images every 64 s, the image triggered interval could be longer than 64 s when the event has also been 
triggered by the rate trigger during the image triggered interval. 
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FIG. 6. — The photon index versus the photon fluence in the 15-50 keV band in the image interval of the 64s image triggered GRBs. The spectral model is a 
simple power-law model. 



is different from their minimum peak flux for the total observed duration, shown in Figure 5. 

We found 33 bursts (out of the 72 image triggered GRBs) with no optical counterpart observed, which are listed in Table 1, by 
looking through the Gamma-ray bursts Coordinates Network (GCN) circulars. These may be candidates for being high-redshift 
bursts with z > 6, and might include Pop. Ill GRBs. Their Tgo durations are relatively long, but not as long as the <~ 1 day, 
predicted for Pop. Ill GRBs in our model. These Tgo durations are just the time intervals during which BAT was able to detect 
90% of the photons from the sources, which is not necessarily the same as the real duration of the bursts, although it provides 
a useful uniform measure approximating this quantity. If the burst flux is marginally above the BAT threshold initially and it 
gradually declines, the Tgo duration could be much shorter than the real duration. 

A bright external shock emission of a Pop. Ill GRB would trigger the BAT even if the prompt emission flux is below the 
detection threshold. We calculated the 64 s photon fluences of the external shock emission in the 15-50 keV band for the cases 
of Figures 2 and 3 and obtained S p h — 0.33 ph cm -2 and ~ 1.5 ph cm -2 , respectively. The photon fluence in the latter case is 
marginally above the effective threshold <~ 1 ph cm -2 . For the case of 1 +z = 10 (Figure 4), we obtained 5 P h — 2.0 ph cm" 2 for the 
negligible pair production case and S p h — 7.5 ph cm" 2 for the significant pair production case. These indicate that BAT will be 
triggered by the external shock emission in some cases. We have a rough relation 5 p h oc e e L^ f{p)dj} (see the end of Section 3.3). 
When the progenitor mass is smaller, leading to L; so three times smaller than the case of Figure 2, the external shock emission for 
the case of e e -\ ~ f(p) ~ 1 and 1 +z > 10 will not trigger BAT (where the pair production is negligible for no < 10 4 ). In the next 
section, we use a specific model of the Pop. Ill GRB prompt emission to examine their detectability. 

4.2. Prompt photospheric emission 

Pop. Ill GRB jets are likely to be dominated by Poynting-flux, as discussed in Section 2. The prompt emission mechanism of 
Poynting-dominated GRB jets (as opposed to the afterglow, on which we concentrated thus far) has been actively discussed in the 
literature (e.g., Thompson 1994; Meszaros & Rees 1997b; Spruit et al. 2001; Lyutikov 2006). The jet may have a subdominant 
thermal energy component of electron-positron pairs and photons, so that the emission from the photosphere can be bright. In 
addition to this, above the photosphere, the magnetic field could be directly converted into radiation via magnetic reconnection 
or the field energy could be converted into particle kinetic energy which can produce non-thermal radiation via shocks. The 
existence of the latter emission components is uncertain and they are currently difficult to model. Thus, for simplicity we focus 
on the photospheric emission, which is essentially unavoidable. Such photospheric emission models of the prompt emission are 
viable also for baryonic jets, which could work for Pop. I/II GRBs (e.g., Meszaros & Rees 2000; Rees & Meszaros 2005; Ioka et 
al. 2007; Toma et al. 2010). MR10 developed the Poynting-dominatedjet model of Meszaros & Rees (1997b) for a Pop. Ill GRB 
jet, and estimated the luminosity and temperature of the photospheric emission. Here we recalculate this emission component, 
taking into account the collimation of the outflow. 

Let us assume for simplicity that the opening angle of the jet is roughly constant from the base of the jet, r; = IghRh — 
9 Ax 10 7 ghMf,,2.5 cm where gh is a numerical factor, out to the external shock region. The isotropic-equivalent luminosity of the 
jet is given by Lj so = L B z(2/# 2 ) — 4.4 x 10 53 (al/a-ip^Mj^.st^Oj 2 ^ erg s" 1 . Denoting by a the ratio of the Poynting energy flux 
and the particle energy flux at the base, the comoving temperature of the flow is estimated as T( = (Li S0 /[(l + a^TrrjcaTj]) 1 ^ ~ 
5.3 x 10 5 [(1 + a)/lQT l l A L l Jl 6 r~ l ^ 2 T~ l 112 eV/Jt, where L 53 .6 = iiso/10 53 ' 6 erg s" 1 , r h% = r,/10 8 cm, and V, is the bulk Lorentz factor 
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TABLE 1 

Prompt emission parameters for 33 image trigger GRBs with no optical counterpart. 790 is 790 duration of the burst. S(15-150) is 

THE FLUENCE IN THE 15-150 KEV BAND. P 1 ^ 15- 150) IS THE 1-S PEAK PHOTON FLUX IN THE 15-150 KEV BAND. F IT ( 15-50) IS THE ENERGY FLUX 

MEASURED IN THE IMAGE TRIGGERED 64 S INTERVAL IN THE 15-50 KEV BAND. PL PHlNDEX IT IS THE PHOTON INDEX BASED ON A SIMPLE 
POWER-LAW MODEL TO THE SPECTRUM OF THE 64 S IMAGE TRIGGERED INTERVAL. 



GRB 


J90 


S(15-150) 


PpJ(15-150) 


F IT (15-50) 


PL PhIndex IT 




(s) 


(erg cm" 2 ) 


(ph cm" 2 s" 1 ) 


(erg cm" 2 s" 1 ) 




0507 14B 


46.9 


(5.3 ± 1.0) x 10"' 


0.6±0.2 


(6.0 ± 1.0) X 10"' 


2.5 ±0.4 


050803 


88.1 


(2.2 ±0.1) x ur" 


1.0±0.1 


(4.2 ± 0.8) x 10"' 


1.6 ±0.3 


050916 


80.0 


(1.1 ±0.2) x ur" 


0.5 ±0.2 


(4.3 ± 0.7) x 10"' 


1.8 ±0.3 


050922B 


156.3 


(2.4 ±0.4) x 10"" 


1.1 ±0.4 


(7.2 ± 1.0) x 10"' 


2.0 ±0.3 


051001 


190.6 


(1.8 ±0.1) x ur 6 


0.5 ±0.1 


(2.6 ± 0.6) x 10"' 


2.0 ±0.3 


051213 


71.1 


(8.o± 1.0) x ur 7 


0.5±0.1 


(3.4 ± 0.6) x 10"' 


1.6 ±0.3 


051221B 


39.9 


(9.1 ± 1.0) x ur 7 


0.6±0.2 


(2.4 ± 0.5) x 10"' 


1.7±0.3 


06021 1A 


118.2 


(1.6 ± o.i) x ur 6 


0.4±0.1 


(3.0 ± 0.3) x 10"' 


1.4 ±0.2 


060319 


8.9 


(2.4 ±0.3) x lO" 7 


1.1 ±0.1 


(2.3 ± 0.4) x 10"' 


2.0 ±0.4 


060413 


117.3 


(3.6±o.i) x ur 6 


0.9 ±0.1 


(2.3 ± 0.5) x 10"' 


1.6±0.3 


060427 


62.0 


(5.0 ±0.9) x ur 7 


0.3±0.1 


(3.1 ±0.5) x 10"' 


2.0 ±0.4 


060516 


161.2 


(l.o ±0.2) x ur 6 


0.4 ±0.2 


(3.4 ± 0.6) x 10"' 


2.0 ±0.3 


060728* 


- 


(2.2 ±0.7) x 10" 7 


0.08 ± 0.02 


(1.2 ±0.4) x 10"' 


1.4 ±0.4 


060923C 


67.4 


(1.6 ±0.2) x ur 6 


0.9 ±0.3 


(1.1 ±0.1) x 10" 8 


2.2 ±0.2 


061027* 


- 


(2.9 ±0.9) x 10" 7 


0.2±0.1 


(2.2 ± 0.7) x 10"' 


1.9 ±0.5 


061028 


105.6 


(9.5 ± 1.7) x lO" 7 


0.6 ±0.2 


(3.0 ± 0.6) x 10"' 


1.9 ±0.4 


070126* 


- 


(1.3 ±0.4) x 10" 7 


0.18 ±0.04 


(1.2 ±0.4) x 10"' 


1.9 ±0.5 


070429A 


168.0 


(9.4 ± 1.4) x 10" 7 


0.4 ±0.1 


(3.9 ± 0.6) x 10"' 


2.2 ±0.3 


070520A 


71.0 


(5.o ± 1.0) x ur 7 


0.3±0.1 


(2.5 ± 0.5) x 10"' 


1.9 ±0.3 


070704 


377.6 


(5.9 ±0.3) x lO" 6 


2.0±0.1 


(1.1 ±0.1) x 10" 8 


2.2 ±0.2 


070920A 


51.5 


(5.2 ±0.7) x 10" 7 


0.3 ±0.1 


(2.9 ± 0.4) x 10"' 


1.6 ±0.2 


071018* 


288.0 


(l.i ±0.2) x ur 6 


< 0.2 


(1.0 ±0.2) x 10"' 


1.6 ±0.3 


071021 


228.7 


(i.4±o.2) x ur" 


0.6±0.1 


(2.6 ± 0.5) x 10"' 


1.8±0.4 


07 1028 A* 


33.0 


(3.3 ±0.6) x HF 7 


< 0.4 


(3.0 ± 0.5) x 10"' 


1.8 ±0.2 


080207* 


> 290 






(6.7 ± 0.5) x 10"' 


1.3 ± 0.1 


080325* 


162.8 


(4.9 ± 0.4) x ur 6 


<2.0 


(1.4 ±0.3) x 10" 8 


1.5 ±0.3 


081017* 


> 320 


(1.4 ±0.2) x ur" 


0.2±0.1 


(1.9 ±0.3) x 10"' 


1.6 ±0.2 


081022 


157.6 


(2.6 ±0.2) x lO" 6 


0.6±0.1 


(7.5 ± 0.6) x 10"' 


1.4 ±0.1 


090308 


25.1 


(2.2 ±0.5) x lO" 7 


0.3±0.1 


(3.1 ±0.5) x 10"' 


2.4 ±0.4 


090401 A 


117.0 


(1.12±0.03) x 10" 5 


11.0±0.4 


(6.3 ± 0.9) x 10"' 


1.6 ±0.2 


090419* 


460.7 


(2.7 ±0.3) x 10"" 




(3.3 ± 0.5) x 10"' 


1.3 ±0.2 


090807A 


146.4 


(2.2 ±0.2) x Hr" 


0.7 ±0.2 


(8.6 ± 0.9) x 10"' 


2.1 ±0.2 


091104 


107.1 


(7.6 ± 1.2) x 10" 7 


0.4±0.1 


(3.0 ± 0.5) x 10"' 


1.7 ±0.3 



t Possible GRB. 

Uncomplete data to measure the fluence and peak flux. 
*Poor statistic to measure the 1-s peak flux form the spectrum. 



of the flow at the base. The dynamics of the flow while it is optically thick is governed by energy conservation (L lso = const.), 

9 /3 

entropy conservation (r z T V = const.), and the MHD condition for the flow velocity to be close to the light speed (the lab-frame 
field strength B oc r" 1 ). The last condition indicates that the Poynting energy is conserved, and so the particle energy is also 
conserved, i.e., r 2 r' 4 T 2 = const, and a = const. Then we have T' oc r -1 and r oc r. At the photosphere radius r = r a where the 
electron-positron pairs recombine, the temperature is given by kT' a ~ 17 keV. This leads to 

r fl ~3.1 x 10 9 [{l+a)/Wr 1/ *L% 4 6 r^r; 1/2 cm, ( 12 ) 

r fl ^3i [(i+ayior^^fr 1 / 2 . (13) 

The observed temperature and the bolometric energy flux of the photospheric emission are then 



F a '. 



1+z 

i26 [(l + CT )/10]- 1 / 4 L^ 4 6 r/";8 /2r / 1/2 [( 1 +^)/ 2 0]" lke V, (14) 
W(l+<7) 



4nd 2 L 



-7.1 x lO" 9 [(l + ( 7)/10]- 1 L 5 3.6fl'Z 2 2oergcm- 2 s" 1 . (15) 
The photospheric photons can be scattered by MHD turbulence or Alfven waves, induced by e.g., the interaction of the jet with the 
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stellar envelope, into a power-law spectrum extending up to comoving photon energies ~ m e c 2 (Thompson 1 994) 10 . The emission 
from the photosphere of Pop. Ill GRB would thus have a black-body peak, and may have a non-thermal tail extending (in the 
observer frame) to photon energies e <~ T a m e c 2 /(l+z) ~ 1 MeV. This photospheric emission may be detected until f b s = f<i, bs- 
We plot this in Figures 2, 3, and 4 with dotted lines for the case of 1 + a = 10 and L53.6 = 0,8 = T/ = 1. 

We calculated the photon fluences in 64 s in the 15-50 keV band for the above parameter sets as S p h ~ 1 .4 ph cm" 2 for the case 
of 1 +z = 20 and ~ 1 .3 ph cm -2 for the case of 1 +z = 10. Thus this emission can be marginally detected by the image trigger of 
BAT. The value of a is highly uncertain, similar to /?» and 6j for a specific value of M* . For L{ so / (1 + a) three times smaller than 
the above case (and similar values of r\ and Ff), the photon fluence is calculated as < 1 ph cm -2 both for 1 +z = 20 and 1 +z = 10, 
and thus BAT is not expected to be triggered. 

4.3. Pop. Ill GRB Rate 

The Pop. Ill GRB rate is largely uncertain, and has only been inferred from theoretical models. The observed rate of Pop. Ill 
GRBs originating between redshifts z and z+dz is computed by 

dN ohs 1 dV 

™ = ^(z)W«)t^?- < 16 > 
dz 1+zdz 

where tjj*(z) is the Pop. Ill star formation rate (SFR) per unit comoving volume, ?7grb(z) is the efficiency of the GRB formation, 
P(z) is the detection efficiency, i.e., the ratio of the Pop. Ill GRBs which would be detected by a specific instrument out of the 
entire number of Pop. Ill GRBs, and dV /dz is the comoving volume element of the observed area per unit redshift. The additional 
factor 1/(1 +z) represents the cosmological time dilation effect. 

The factors ip*(z) and ?7grb(z) are both highly uncertain for the Pop. Ill VMSs as well as for the Pop. I/II stars. Just for a concrete 
discussion, here we use i[>*(z) predicted by using the extended Press-Schechter formalism and the current observational results on 
the SFR for z < 6 (Bromm & Loeb 2006) (see also Naoz & Bromberg 2007). For Pop. I/II stars, they assumed a z-independent 
7?grb ~ 2 x 10" 9 Mq 1 to set the total Pop. I/II GRB rate observed by Swift BAT to be ~ 90 yr" 1 . (Note that they assumed that 
Swift BAT covers Air of the sky, while it actually covers only ~ 27r/3 of the sky. Thus the GRB formation efficiency should be 
normalized as ~ 6 times their adopted value, 7?grb — 1.2 x 10~ 8 Mq.) Their result predicts the Pop. I/II GRB rate observed by 
BAT at z > 6 to be ~ 10 yr" 1 . Since a small fraction ~ 25% of GRBs detected by BAT have redshifts determined, because of bad 
conditions for optical and near-IR observations (cf. Fynbo et al. 2009, see also our implication from the image triggered GRBs in 
Section 4. 1 and in Table 1), the predicted rate of Pop. I/II GRBs with z determined would be ~ 3 yr" 1 . This is somewhat higher 
than the current observed rate, 0.6 yr" 1 , i.e., 3 GRBs (GRB 050904, GRB 080913, and GRB 090423) during the 5-yr operation of 
Swift, but their model of V>*(z) and 7/grb is not interpreted as unacceptable, taking into account the uncertainties of the theoretical 
calculations and the poor statistics of the current observed data. For Pop. Ill stars, they assumed the same z-independent i]grb 
as Pop. I/II stars and computed the nominal Pop. Ill GRB rate observed by BAT to be /?bat ~ zdN^ s /dz ~ 0.03 yr" 1 for bursts 
around z ~ 20 and <~ 0.3 yr" 1 for bursts around z <~ 10. 

The detection efficiency P(z) is computed by assuming the GRB luminosity function and the detection threshold. Bromm & 
Loeb (2006) assumed that the Pop. Ill GRB luminosity function is the same as that of Pop. I/II GRBs and took a BAT detection 
threshold of / p h,iim — 0.2 ph cm" 2 s" 1 . We can assume, however, that the Pop. Ill GRBs have a different luminosity function, with 
a brighter membership than Pop. I/II GRBs, and take the effective BAT threshold deduced by our analysis of the image trigger 
bursts (see Section 4. 1), in which case Rbat can be larger than the above estimate. 

In order to estimate /? B at in our Pop. Ill GRB model, let us first presume the observed rate (for <~ 27r/3 of the sky) if we would 
detect the entirety of the Pop. Ill GRBs originating around z ~ 10-20 without considering the detection efficiency, R w , based on 
the calculations of Bromm & Loeb (2006). The detection efficiency is computed by 



P(z)= / P (L ph )dL ph , (17) 

where L p h is the isotropic-equivalent photon luminosity of a burst, and p(L ? h) is the luminosity function normalized as 
/ °° p(Lp h )dLp h = 1. They adopted L p h,ii m = 47rJ|/ p h,iim - lx 10 60 ph s" 1 for bursts around z ~ 20 and ~ 2 x 10 59 ph s" 1 for 
bursts around z ~ 10. The detection efficiency is then roughly estimated to be ~ p(L v h.Um)Lph.Um/ p(Lp)L p ~ 0.2 for z ~ 20 and 
<~ 0.5 for z ~ 10, where L p provides the peak of the function of /?(L p h)L p h. Thus, we obtain R w <~ R BA t/0.2 <~ 0.2 yr" 1 for z ~ 20 
and - /?bat/0.5 - 0.6 yr" 1 forz ~ 10. 

According to our study in Sections 4. 1 and 4.2, the photospheric prompt emission and/or the external shock emission of Pop. Ill 
GRBs can be detected by the BAT image trigger for our fiducial set of parameters, Md.2.5 — ^,2.5 —R*,\2 — — — fip) — 
(l + cr)/10~ 1 and l+z> 10, while not detected for L; so just three times smaller than that for the above set of parameters. Thus, 
if the distribution of the parameters of the Pop. Ill VMSs clusters around our fiducial set of parameters, BAT would detect, say, 
about half of the Pop. Ill GRBs out of the whole Pop. Ill GRBs, i.e., the detection rate could be /? B at ~ ^w/2 ~ 0.1 yr" 1 for 
z ~ 20 and <~ 0.3 yr" 1 for z ~ 10. These values have large uncertainties, but imply that the 5-yr operation of Swift so far may 
already have detected a Pop. Ill GRB, or may detect it in the near future, if the factor ip*(z) x t?grb is given as above. 

10 If the power-law spectrum extends to energies much higher than ~ m e c 2 , e.g., due to magnetic dissipation, as argued in MR10, copious pair formation 
would ensue, which would form a new (pair) photosphere at a larger radius (e.g., Rees & Meszaros 2005). 
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5. LATE-TIME RADIO AFTERGLOWS 

We have focused so far on the high-energy emission just before and near the beginning of the external shock self-similar phase 
(f b s ~ td,obs, see Figure 1) to constrain the physical parameters £i S0 and n of Pop. Ill GRBs, and to examine their detectability by 
BAT. In this section, we argue that the radio afterglows of Pop. Ill GRBs in the self-similar phase (f 0DS > frf )0 bs) can be so bright 
that they also provide powerful tools for constraining the event rate. At ?d,obs> the flux of the external shock emission shown in 

Figure 2 at v a ~ 230E^ 6 6 t~^ 2 n l Q /6 [(l +y)/3.7]" 1 / 3 [(l +z)/20]~ 1 GHz is given by 

F Va ~ 2.6 e^ 4 4 /3 6"o I/12 [(l +i0/3.7]- 5 / 6 [(l +z)/20]d£ a Jy. (18) 

-1/2 

The break frequency v a decreases (as f obs ' ) but the flux at v a stays constant after fd,obs> at least until the epoch when the jet 
effects become significant. This indicates that the Pop. Ill GRB afterglows can be very bright radio sources, despite their large 
distances. 

We briefly compute the light curves at various frequencies in the radio bands, 100 GHz, 5 GHz, 1 GHz, and 70 MHz, in our 
fiducial case, E il fl = t d ^ = n = £b-2 = = f(p) = = 1- As discussed above, the temporal evolution of the characteristic 
quantities of the external shock in the self-similar expansion phase can be obtained by replacing t d by the variable f obs and 
taking the other parameters as constant in the equations of the general afterglow model (shown in Appendix). The jet effects are 
significant when the Lorentz factor of the shocked fluid is T ~ 8~ l , i.e., at the observer's time 

fe.obs - f A ™ m J - 5 J d+z)-9.9x 10 2 E^n-^ef, [(1 + z )/20] day. (19) 

After tg t0 b s , the temporal evolution of the characteristic quantities is obtained by replacing £i S0 by t ~ b ' s and td by f ODS , respectively, 
and taking the other parameters as constant (Sari et al. 1999). Such an evolution is derived under the assumption that the shocked 
fluid expands sideways rapidly after fg.obs- Recent detailed hydrodynamic simulations (Zhang & MacFadyen 2009; Granot 2007) 
have shown that the shocked fluid only undergoes a slow sideways expansion, while the afterglow light curves can be still 
approximated by those predicted by Sari et al. (1999) especially in the radio band. (In the optical and X-ray bands, the spectrum 
is so soft that the limb-brightening effect is significant. Then the brightest portion at an angle 9 <~ T~ l from the line of sight 
becomes missing at f^,obs> which causes a steeper light curve than f"^.) 

In an analytical treatment (e.g., Sari et al. 1999) the rapid sideways expansion of the fluid would lead to a nearly spherical 
(Oj <~ 1) configuration of the shocked fluid leading to the non-relativistic Sedov-von Neumann-Taylor (SNT) solution, starting 

around observer's time ? S NT,obs - [E^O] /(Aimm p c 5 )] l ^{\ + z) ~ 9.9 x 10 4 E^n^^ef^l +Z V 2 °] da y- However, the detailed 
hydrodynamic simulations (Zhang & MacFadyen 2009) show that the much slower sideways expansion results in the later start 
of the SNT phase. The starting time is shown to be a few times ?NR,obs 



1/3 

\Aimm p c 5 



( — 3 (l + z)^4.6x 10 5 £' /3 .«- 1/3 



Ws * . 5 (1 + z) 4.6 x 10 s E^>o ' [(! +^)/ 2 °] da y> (2°) 



where the shocked fluid is still highly collimated. The temporal evolution of the characteristic quantities in the SNT phase is 

obtained by using r oc f 2 ^ 5 and v oc t Q ^ 5 , where v is the velocity of the shocked fluid, and the internal energy density is oc v 2 . This 
results in a light curve with a shallower decay at v > v a (or a steeper rise at v < v a ) in the SNT phase, compared to the decay 
in the prior phase. Therefore, the assumption of the slow sideways expansion and the late SNT phase based on the numerical 
simulations leads to smaller number of off-axis observers (due to the high collimation) and dimmer radio fluxes than that of the 
rapid sideways expansion based on the approximate analytical arguments. We here take the former, conservative assumption of 
the slow sideways expansion. 

For calculating the radio emission we need to compute only the synchrotron emission of the original electrons. At late times, 
the number of pairs is typically small, and the pair emission is negligible. The SSC emission is not relevant in the radio bands. 
The results are plotted in Figure 7. Here we confirmed the assumptions that all the electrons in the shocked region remain 
relativistic at least until f D b s = 10 6 day. 

Figure 7 shows that the radio afterglows of Pop. Ill GRBs can be very bright with a very long duration. These could have been 
detected as quasi-steady point sources by the radio survey observations. As far as we know, the current largest radio survey data 
is based on the Very Large Array (VLA) FIRST survey (White et al. 1997), which observed a large area mainly around the north 
Galactic cap at 1 .4 GHz, covering ~ 1 /5 of all the sky. 1 1 For the threshold, <~ 6 mJy at 1 .4 GHz, a Pop. Ill GRB radio afterglow at 
1 +z = 20 (at 1 +z = 10) can be observed for f r d i0 b s ~ 300 yr (for f r j. b s ~ 200 yr). Therefore, if we denote by R w .4n the all-sky Pop. Ill 
GRB rate in unit of yr" 1 , the number of the Pop. Ill radio afterglows that would have been detected in that survey is estimated to 
be - {R wM /5)t rd , ohs ~ 60 (R wM /l yr" 1 )(f r < J , o bs/300 yr) for bursts at 1 +z ~ 20, and - 200 (R wM /4 yr 1 )(f^ !Obs /200 yr) for bursts 
at 1 +z ~ 10. A detailed analysis of the FIRST data would thus provide a powerful constraint on R w ^ of Pop. Ill GRBs such 
as discussed in this paper (even no sources like our model calculations would provide an upper limit on the rate). The Pop. Ill 
radio sources could have X-ray counterparts. In our model the late time X-ray afterglow is dominated by the SSC component, 

11 Levinson et al. (2002) and Gal- Yam et al. (2006) did not find any radio transient sources like GRB afterglows with timescales of significant flux changes 
~ 5 yr by comparison between the NVSS (spanned over 1993-1996) and FIRST (1994-2001) catalogs, which effectively cover ~ 1/17 of the sky. This indicates 
thatR^ < 17/5 ~ Syr- 1 . 
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FIG. 7. — Radio light curves at frequencies, 100 GHz (dot-dashed line), 5 GHz (dotted line), 1 GHz (solid line), and 70 MHz (dashed line), of a Pop. Ill GRB 
at 1 +z = 20 with typical parameters £57.6 = ?rf.4 = «o = e B,-2 = £e,-l = f(p) = 0j,-l = 1 ■ The important observer's times are the jet duration td, bs — 2.3 day, the jet 
break time tg obs ~ 9.9 X 1() 2 day, and the time when the shock becomes non-relativistic ?NR,obs — 4.6 X 10 s day. 



whose flux in the 0.3- 10 keV range can be > 10" 15 erg cm" 2 s" 1 for f obs < 100 yr (which can be detected by Chandra X-ray 
Observatory). 

The above predicted numbers of the Pop. Ill GRB radio afterglows are just based on the values of f ri i, bs calculated for our 
fiducial set of parameters. In order to obtain more realistic numbers taking into account the distributions of the model parameters, 
we require to calculate the radio light curves for large ranges of parameter values. However, the radio light curves may highly 
depend on the parameters £i S0 ,^,n,efl,e e , and f(p). In other words, there are many patterns of light curves depending on the 
orders of characteristic times t a=c ,t m=a , and t m=c , which are the times when v a crosses v c , when v m crosses v a , and v m crosses 
v c , respectively, as well as td,te, and ?nr (we have tj < t$ < t m=a < t m=c < t a=c < /nr for our adopted parameter set). A thorough 
model analysis with different parameter sets would deserve another future work. 

Recent radio transient searches without primary detections at any other frequencies have begun to open new observational 
frontiers (e.g., Niinuma et al. 2007; Bower et al. 2007; Ofek et al. 2010), but currently they do not appear to have good potentials 
for detecting Pop. Ill GRB radio afterglows. The interferometric drift-scanning observation with the Waseda Nasu Pulsar Obser- 
vatory in Japan may scan a large area at 1.4 GHz (Niinuma et al. 2007), but its detection thereshold, <~ 0.3 Jy, is not sufficiently 
high for searching our typical Pop. Ill radio afterglows. Bower et al. (2007) reported results of a survey for transient sources by 
using archival data obtained with VLA observation of a single field at 5 or 8.4 GHz, but the field is very small, < 0.1 deg 2 . 

The absorption line at 21 cm (associated with transitions between the hyperfine levels of the hydrogen atom) seen in the 
continuum radio spectrum of high-redshift sources would be one of the promising tools to reveal the reionization histroy of 
the universe (Tozzi et al. 2000; Shaver et al. 1999). GRB radio afterglows have been considered as candidate backlighting 
sources, but the bursts with the usual or slightly higher total energies have been found to be too dim for this aim (Furlanetto 
& Loeb 2002; Ioka & Meszaros 2005). However, the high-redshift GRBs arising from Pop. Ill VMSs, such as discussed in 
this paper, can emit a sufficiently bright radio afterglow to be of interest. For our fiducial parameter values, the radio flux at 
v = 1420 MHz/(l +z) is estimated to be ~ 2.4 mJy (with peak time f p o b s — 850 yr) for 1 +z = 20 and ~ 6.4 mJy (with peak time 
f^.obs ~ 430 yr) for l+z = 10, which are comparable to the lower bound on the radio flux for detection of the 21 cm absorption 
line, ~ 1 — 10 mJy (Ioka & Meszaros 2005). Thus the 21 cm absorption line could be marginally detected in the Pop. Ill GRB 
radio afterglows. A radio survey with telescopes like Low Frequency Array (LOFAR) 12 could detect such radio afterglows 
and determine their redshifts by their 21 cm lines themselves. The predicted number of detections is <~ /? w ,47r(fi/4vT)? J , )0 bs ~ 
1 (tfwWl yr 1 )(Q/50 deg 2 )(^, obs /850 yr) for bursts at 1 +z ~ 20, and - 1 (R wM /4 yr" 1 )(O/20 deg 2 )(f p>obs /430 yr) for bursts at 
1 +z ~ 10, where Q is the solid angle of the survey area. Note that the jets will keep collimated even around ~ fNR.obs due to the 
slow sideways expansions (Zhang & MacFadyen 2009), so that off-axis viewings of the afterglows are too dim to detect. 13 

6. SUMMARY 

12 http://www.lofar.org. 

13 If the jets were instead to undergo a rapid sideways expansion, as suggested by the approximate analytical treatment of the transition to the non-relativistic 
regime (e.g., Sari et al. 1999), the peak fluxes and peak times are slightly larger than the above estimates, i.e., F p ~ 4.3 mJy with f^ obs — 1300 yr at 1 +z = 20, 
and F p ~ 1 1 mJy with t p obs ~ 630 yr at 1 +z = 10. Since t p (, s > (sNT.obs an d thus the fluids would be nearly spherical at the peak time, the off-axis observers 
can detect the emission. Therefore the predicted number of detections would be much larger, by a factor of 0y 2 = I0 2 d~j 2 _ 1 . 



14 



Toma et al. 



Pop. Ill GRBs may have an isotropic-equivalent energy Ei so > 10 erg and a cosmological-rest-frame duration td > 10 s. We 
have calculated the external shock emission spectrum at f^obs in detail based on the standard model (Meszaros & Rees 1997a; Sari 
et al. 1998; Sari & Esin 2001; Nakar et al. 2009). This model can explain many of the X-ray/optical/radio afterglows detected so 
far (e.g., Panaitescu & Kumar 2002; Liang et al. 2007) as well as the very early high-energy afterglows detected by Fermi LAT 
(e.g., Kumar & Barniol Duran 2009; De Pasquale et al. 2010; Corsi et al. 2010). 

We found that the external shock emission at td, bs can be detected by Fermi LAT and Swift XRT/BAT, whose flux leads to a 
constraint on E- 1%0 by using the source redshift (and distance) that will be determined by the observation of the Lya drop-off in 
the IR band. The detection of a burst at z > 10 with E\ so > 10 57 erg and t d > 10 4 s would be a very strong indication that this is a 
GRB arising from a Pop. Ill VMS with M* > 300 Mq. This indication should be complemented with the constraint on the metal 
abundances in the surrounding medium through high resolution IR and X-ray spectroscopy. 

One of the important findings of the present study is that the 77 self-absorption break at energy e 77 in some cases of the 
external shock emission spectrum of a Pop. Ill GRB is expected to be observable in the LAT energy range. Given that the prompt 
emission at f^,obs does not hide the external shock emission in the LAT range and e~ n is well below the values of the EBL 77 
cutoff energy expected for practically all EBL models with the determined source redshift, we have shown that the flux and the 
energy e 77 , together with Swift XRT data, can lead to a constraint on the value of the external density n. The constraint on n 
would provide invaluable information about the environment and the radiative feedback processes of Pop. Ill stars. 

Putting constraints on n from the 77 self-absorption break is a fairly new method, while constraints on n from the multi- 
wavelength observations of the late-time afterglows of Pop. I GRBs are common (e.g., Panaitescu & Kumar 2002; Wijers & 
Galama 1999). Here we have discussed the conditions under which these methods can be used in Pop. Ill GRBs. One caveat is 
that we have assumed that the electron acceleration processes works uniformly in the emitting region for calculating the afterglow 
spectra. As discussed below Equation (A19) in Appendix, it is also possible that the electron acceleration works only near the 
shock front, where e 77 may be far above the LAT energy range, and Sm depends on the unknown upstream magnetic field strength, 
which do not allow us to constrain n. 

In order to identify the direction to GRB afterglows on the sky for observations with Swift, Fermi, and the IR telescopes, the 
emission needs to be high enough to trigger the large field instrument BAT. A reasonable trigger threshold of Swift BAT can be 
estimated by focusing on the image trigger mode, since it is based on a criterion with a fixed time-scale and energy band, while 
more general BAT threshold including the usual rate trigger is too complicated to estimate. The image trigger mode may be 
suitable for detecting weak and less-variable bursts like very-high-redshift bursts. We have used the samples in the BAT2 catalog 
and deduced the detection threshold of the BAT image trigger to be ~ 1 ph cm" 2 for the 64 s interval in the 15-50keV band. 
We have calculated the prompt photospheric emission flux of Pop. Ill GRBs with typical parameters, and shown that these can 
be marginally detected by BAT. The external shock emission can also trigger BAT without the prompt emission trigger in some 
cases. 

We have also briefly shown that the Pop. Ill GRB late-time radio afterglows can be very bright. For our fiducial parameters, 
the peak flux at 1 GHz is ~ 140 mJy, which could be identified in the VLA FIRST survey data. This survey covered a sufficiently 
large area, which would provide a powerful upper bound on the rate of the Pop. Ill GRBs. The peak flux of the late-time 70 MHz 
radio afterglow for our fiducial parameters is ~ 2.4 mJy, in which 21 cm absorption lines could be detected. This would provide 
a measure of the neutral hydrogen fraction in the IGM around the Pop. Ill star. 

Putting constraints on the properties of Pop. Ill stars has recently become of great importance in modern cosmology. Planned 
IR surveys will be able to probe Pop. Ill stars. However, it is difficult to distinguish between a single Pop. Ill VMS and a cluster 
of less massive Pop. Ill stars. Thus, the detection of GRBs with very high E\ so and very long td could provide critical, 'smoking 
gun' evidence for the existence of VMSs. Multi-wavelength observations of such GRBs with Swift, Fermi, and ground-based IR 
and radio telescopes should provide us with invaluable information on Pop. Ill stars and their environments. 

We thank D. N. Burrows, A. D. Falcone, D. B. Fox, A. Gal- Yam, S. Gao, K. Murase, S. Naoz, and the anonymous referee for 
useful comments. We acknowledge NASA NNX09AT72G, NASA NNX08AL40G, and NSF PHY-0757155 for partial support. 
PM is grateful for the hospitality of Fermilab and the Institute of Astronomy, Cambridge University, during part of this project. 



The afterglow emission spectrum at the time td is determined by the radius and the Lorentz factor Yd of the shocked fluid 
at this time, which are given via the relations £i S0 ~ 4nr d T d nm p c 2 and ~ cY 2 d td- Here we have assumed for simplicity that the 
circumburst medium density is uniform, n = 1 «o cm -3 . These two equations lead to 
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where £57.6 = £i S o/10 57 6 erg and td,4 = f^/10 4 s. The magnetic field strength in the shocked region scales as 

B ~ {^€ B nm p c 2 ) l l 2 T d ~ 3.8 eJ^E^t^J*^ G > 



(A3) 
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where eg-2 = eg/lO -2 . The minimum injected electron Lorentz factor is 

m e p-1 

where e e _j = e e /10 _1 and /(p) = (13/3)(p — 2)/(p — 1). We have assumed that all the electrons are accelerated to a power-law 
spectrum dn/dj e oc 7j p for j e >j m . 

The accelerated electrons will lead to synchrotron and SSC emission. The radiative cooling timescale in the comoving frame 
is t' c = 6irm e c/[<jT7eB 2 (l + T( 7e ))], where T(7 e ) is the luminosity ratio of the SSC to synchrotron emission for electrons with 7 C , 
while the comoving dynamical timescale is t' d = T d t d . Thus the electron Lorentz factor above which the radiative cooling is more 
significant than the adiabatic cooling is 

6nm e c j _j _3/8 i/8 -5/8 

All the injected electrons are radiatively cooled within the dynamical timescale if j c < 7m , that is if 

e B ,-2nJ\\+Y{ lc )\ > 1.4 x 10" 2 e^E^J/OO -1 - (A6) 
This condition can be rewritten using Eqs. (3) and (6) as 

^ d %K^Kn n f lU a ^' 2 ^-^ 1 ' 2 > lO-^e-U/Xpr'fl +r(7c)] _1 . (A7) 

Since the masses and M/, are expected to be positively correlated with the radius R*, this condition is found to be satisfied for 
reasonable parameter ranges. We thus focus on the fast-cooling regime, -f c < -f m for the emission at t d . 
The peak energies of the eF £ spectra of the synchrotron and SSC emission are given by 

£m ~ ^c^JVz ~ 5A <-^ E ^ttfiP)\iy + z)/20]- 1 eV, (A8) 
^ c ~2 7 2 £ m ~ 1.8 x 10 2 e*_ 1 e^ 2 ^2 / V /4 /0') 4 [(l+z)/20r 1 MeV, (A9) 

respectively. If < T d j m m e c 2 /(l +z), the Klein-Nishina (KN) suppression of the SSC emission is not significant (e.g., Nakar 
et al. 2009). This condition is rewritten as eg-2 < 3.2 x 10 3 e~l 1 £ , 5 7 1 6 fj 4 /(/7)~ 6 , which is satisfied for reasonable parameter ranges. 
In the case of negligible KN effects, Y does not depend on -f e for j e < 7m , which is calculated from 

r d [°° j_, dn 2 _ _ p-l _ e e 1 



TrfA. dj e p-2 e B l+Y 



1c 

where r = cr T r d n is the optical depth for the electron scattering. We obtain Y w y/e e /es ~ 3 t^t x €^f\ for the case of e e > e B , 
while otherwise Y rts e e /e B < 1. The SSC spectrum above e^ c has a softening break at 

£KN * 7 g^- * 1.2e7 2 1 e^ 2 2 £ 5 t^o 1/4 /(pr 2 [(l+z)/20]- 1 TeV. (All) 
The fluxes at £ m and at ef^ are given by 



Audi P-H+Y 

~ 1.6 x 10" 9 (l + r)- 1 e e ,_i£ S 7.^j4/(p)d^o er 8 cm " 2 s ~'> (A12) 
4 c F £ s „ c ~T£ mJ F e „,, (A13) 



x 10 29 cm. In the case of significant KN effects, 



where the luminosity distance d L is normalized by the value for 1 +z = 20, 6.7 

i.e., > T d j m m e c /(I + z), the SSC emission is not important and Y is smaller than the above value (see Nakar et al. 2009, for 
details). 

The maximum energy of the electrons is determined by equating the acceleration timescale and the radiative cooling 
timescale. If the electron acceleration occurs in the whole region, we may estimate the (comoving) acceleration timescale as 
t' a = g2irj e m e c I (eB) where g is a numerical factor. For electrons with maximum Lorentz factor -j M , the KN effect is typically 
significant, and T( 7 m) < 1 ■ Then we have ^ M — [3e /(ga T B)] '/ 2 , and its characteristic synchrotron energy is 

3heB 2 T d 9he 2 T d 



£ M— ~. 1m 



4-nm e c l+z 4-irm e c<jTg 1+Z 
~ 1.8 x 10 2 g- l El^f d fn- Xl \{\ +z)/20]" 1 MeV. (A14) 

The energy e a below which the synchrotron self-absorption effect is significant can be estimated by equating the synchrotron 
flux to the blackbody flux of the characteristic electrons in the shocked region, F 6a = [(1 + z) 3 / l d 2 ]2nm e 'y c h(£ a / h) 2 (r^/I ' d ) (e.g., 
Meszaros & Rees 1997a). The Lorentz factor 7^ of the characteristic electrons is given by 7^ whose synchrotron energy is 
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e a in the case of 7a > j c , and otherwise by 7c . Although more detailed calculations can be done by using the self-absorption 
coefficients (e.g., Melrose 1980; Matsumiya & Ioka 2003; Toma et al. 2008), the above approximate derivation is sufficient for 
our current study. Thus we obtain 

j^Yr^E-^^, for 7e < 7fl < 7m , 
la ~ \ 81 (1 + Y)^y 5 iX: »i 7/4l \ ^ la < lc < 7„ ( } 

The characteristic synchrotron energy of electrons with this Lorentz factor is given by 

f 1.5 x lO" 3 (l+y)-V 3 ^/ 6 6f - 1 / 2 „y 6 [(l+ z )/20]- 1 eV, for 7e < 7fl < 7m , 

£a ~ 1 i 1 1A-3/1 , V \ 6/5 e.7/10 -1/2 H/I0 r/1 , . /om -l f . (A16) 

[ 2.1 x 10 (1 +T)e i /_ 2 £ 5 ^ 6 f rf 4 ' n [(1 +z)/20] eV for 7fl < 7c < 7m . 

The high-energy absorption turnover energy e 77 due to the e + e~ pair creation can be estimated as follows (e.g., Zhang & 
Meszaros 2001; Lithwick & Sari 2001). The photons are assumed to be roughly uniform over the emitting region and isotropic in 
its comoving frame, so that e at ~ T 2 d m 2 e c 4 /[e 11 {\ + z) 2 ] is the main target photon energy for the photons with energy £ 77 . Taking 
the spectral luminosity around e = e a , as L £ = L £al (e/e at )~ x , the total number of photons with e > e al is N >£at ~ td J^(L £ /e)de ~ 
L £ jd- By using L £at = \-nd\F £al j{\ +z), we obtain the opacity as r 77 (e 77 ) ~ (<t t /lO)N >Eal /(4Trr d ) = 1, leading to 

e al -9.0 x 10 2 (170)^(1 +Yyhl%l~} 2 E 1 5 %Q + Klf(p) 2 -h (A17) 
where we have taken A = p/2 since e at is typically above e,„. The 77 self-absorption energy is then given by 

e 77 ~ 6.9 (170)^(1 ^^^JE^^n^fipT 2 ^ [(1 +z)/20]" 1 GeV. (A18) 
The condition for a significant absorption, £ 77 < £^ c , is rewritten as 



l-i i _2 £23 



e B ; 2 <(l+T)-? > 38 (\10)^e e _{E- 5 ) 6 tJ'f(p)-^T,. (A19) 

If this condition is satisfied and e e > e B , the emission from the created pairs will significantly affect the observed spectrum. 

Above we have assumed that the photon field is roughly uniform in the emitting region. This may not be valid, however, 
when electrons are in the fast-cooling regime (as is the case for our typical parameters). If the acceleration process of electrons 
works only near the shock front, the emission at e > e m is only produced in a thin layer from the shock front with a width of 
~ (jc/jm)fd/^d in the comoving frame. The annihilation process of the high-energy photons then mainly occurs outside the 
shocked region. In such a case the angles of the directions of motion of two given annihilating photons are typically very small, 
which significantly reduces the pair creation optical depth (Granot et al. 2008). The 77 self-absorption energy e 77 could then be 
much larger than the above estimate by Equation (A18). In addition to this, the maximum synchrotron energy Em is determined 
by the magnetic field in the upstream region, instead of the emitting region (Li & Waxman 2006), which may be smaller than the 
estimate by Equation (A14) by a factor of B u Td/B, where B u is the upstream field strength measured in its own frame. 

An issue which has been a long-standing concern with the model is the possibility that the magnetic field amplification and the 
electron acceleration microphysics may be confined to a region near the shock front. Such small-scale magnetic field will decay 
within a couple of ion skin depths (e.g., Gruzinov 2001; Kato 2005), which may not explain the observed bright afterglows. On 
the other hand, the large-scale magnetic field amplified by macroscopic turbulence created by shock could survive over the whole 
emitting region, and the observed afterglows may be attributed to such large-scale field (Sironi & Goodman 2007). In this case 
the main acceleration process of electrons, e.g., second-order Fermi acceleration, could work uniformly in the shocked region, 
where the photons can be assumed to be roughly uniform and isotropic in the emitting region. If this is the case, e 77 and Sm 
are given by the uniform photon field assumption, Equations (A18) and (A14), respectively. We take this uniform photon field 
assumption in the rest of Appendix and in the main text. 

Case of Negligible Pair Production 

The parameter set of £57.6 = ^4 = «o = £b,-2 = = fip) = 1 satisfies the negligible pair production condition e^ c < e 77 . The 
radius and Lorentz factor of the shocked fluid are given by rj ~ 2.8 x 10 18 cm and ~ 97, respectively. We have Y ~ 2.7, and 
then the characteristic electron Lorentz factors are 7m ~ 4.1 x 10 3 , 7fl ~ 55, and 7c ~ 15. The characteristic photon energies are 
e m ~ 5.4 eV, e a ~ 9.7 x 10" 4 eV, ~ 1.8 x 10 2 MeV (< T d -/ m m e c 2 /(l +z) ~ 10 GeV), e KN ^ 1.2 TeV, e at ~ 2.9 x 10 2 eV, and 
e 77 ~ 21 GeV. The flux normalization is given by e m F £m ~ 4.4 x 10" 10 erg cm" 2 s" 1 , and s^F £ ^ ~ 1 .2 x 10" 9 erg cm" 2 s" 1 . 

The overall spectrum for this case is shown in Figure 2. The synchrotron spectrum is approximately 

1/2 

-M , for £ < e a , 



eF £ ~ e m F £m x { 



for e a <£<£,„, (A20) 

/ ™\"i +1 
\t) , for e m < £ < e M . 
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The SSC spectrum is approximately 



eF E 



(#) (t) ' for £<£fl ' 
, for e a < e < 4°, 

' f ° r £m < £ < £kn ' 



(A21) 



, for e>e KN , 



where we have neglected the segment of e,F e oc e p+1 at e > £ K n just for simplicity (see Nakar et al. 2009, for more details). The 
SSC emission is absorbed by the pair creation within the emitting region at e > e 77 . The spectral shape is given by the intrinsic 
one multiplied by 1 / T 77 (e) oc e~ p l 2 . 



The parameter set of £57.6 = td,A = es-i 



Case of Significant Pair Production 
■ f(p) = 1, no = 10 2 , and e e -\ = 2 satisfies the significant pair production condition 



e 77 < e^ c . The radius and Lorentz factor of the shocked fluid at td are r d ~ 8.9 x 10 17 cm and ~ 55, respectively. The 
characteristic Lorentz factors of the electrons are j m ~ 4.6 x 10 3 , 7 fl ~ 42, and 7 C ~ 0.55, where we have used the value of Y 
calculated below. The cooling Lorentz factor -f c < 1 would imply that the electron energy distribution has a bump at 7^ ~ 1 . 
The synchrotron self-absorption effect may be so strong that the electrons at -f e ~ 1 can be reheated to higher energies. In any 
case, however, the electron energy distribution above j a is not affected, and the synchrotron spectrum below e a is not relevant. 
The characteristic photon energies are e m ~ 21 eV, e a ~ 1.8 x 10" 3 eV, e^ c ~ 9.2 x 10 2 MeV (< r d 7 m m e c 2 /(l +z) ~ 6.5 GeV), 
£kn — 91 GeV, e at ~3.2x 10 3 eV, ande 77 ~6.1 x 1 2 MeV. The flux normalizations are given by e m F €m ~5.8 x 10" 10 ergcm" 2 s" 1 
and e^Fj^ ~ 2.7 x 10" 9 erg cm -2 s _1 . 

The pair creation opacity is larger than unity for photons with ejjjf, so most of the SSC radiation luminosity will be converted 
into e + e~ pairs. If the compactness parameter V in the shocked region is as high as <~ 10 2 , the radiation will be thermalized by 
Compton scattering on the pairs, resulting in a Wien spectrum (Pe'er & Waxman 2004). In our case, 

= <? T e e E /^x 3/ B 1/8 / n N V» 

8irm e c 3 r d r d \Q.2J 576 d ' 4 VlO 2 / ' 



which is high compared to typical cases of Pop. I/II afterglows, but still not as high as 10 2 . We can approximately derive 
non-thermal synchrotron and inverse Compton emission spectra by considering the following cascade process. 

The Lorentz factors of the created pairs are determined by j±m e c 2 w (1 + z)eY~ d l /2. The indices of the injected energy distri- 
bution of the pairs are the same as those of the photons with e > e 77 , being dn±/d^± oc ~f~± 2 for j± tin < 7± < 7±j» oc ^^f^ 1 
f° r l±,p < 7± < 7±,*> an d oc 7/ ^ 2 for y± > 7±,jt, where 



7±,m : 



-77 



(1+z) 



-SC 



2r jm e c 2 ' 2r rf m c c 2 ' ,:c '' i 2r</m e c 2 

For the parameter values adopted here we have 7± m ~ 220, j± jP ~ 330, and 7± ) * ~ 3.2 x 10 4 . The corresponding synchrotron 
energies are e± fn ~ 4.7 x 10~ 2 eV, £± iP ~ 0. 1 1 eV, and £± ^ ~ 1 .0 keV, respectively. Since the synchrotron and SSC cooling is 
very fast, the averaged energy distribution of the pairs over the dynamical timescale is given by 

'7± 2 , for 7±, c <7± <7±,f 

dn± 



(1+z) 



1±,k ^ 



£kn(1+z) 



(A23) 



oc < 



,m ; 

7 ± '7 for 7±, m <7± <7±, P , 

-(p/2)-2 



-5/2 



7± 



(A24) 



-p-(3/2) 



for 7±, P < 7± < 7±,*> 
^7 ± ' - ' for 7 ± > 7 ±jjt . 

The shapes of the synchrotron and inverse Compton emission spectra of the pairs are obtained straightforwardly. The SSC 
emission of the pairs peaks at E*± p ~ 27^ p e±. p ~ 23 keV, and the pairs IC-scatter the original synchrotron emission, which 



peaks at e± 



27j_ p £ m ~ 4.5 MeV. The KN effect is not significant since e 



sc ic 



< T dl ±, p m e c 2 /{\ +z) ^ 4.6 x 10 2 MeV. 



The original electrons IC-scatter the pair synchrotron emission, which peaks at ~ 2j% l £±,p = e± p . 
In order to determine the flux normalization of the various emission components, we define 



t »i e m 


T sc , 


ic f ic 

m e m y 




£mFe m 






AC pIC 







IC £ ±,P^e±j 



£mF E , 



= X, 



AC 



(A25) 



(A26) 



Here we cannot use Equations (A10) and (A13), using instead the following relations for the scattering by the original electrons, 



^sc , 



ylC 

X 



1 



>-2 e B l+y' 



(A27) 
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where Y = Y SC + Y 1C . 
The total number of the pairs can be estimated as 



N± ~ 2t, 



J e 



ix3 rSC 



-de ~ 2L? C t d = 8irdM c 



rSC 



-1/2 



td 
1+Z 



el c (l+z) \e 



-77 

rSC 



-1/2 



SC 



/?-2 y 
p-i i+y 



On the other hand, the number of the original electrons is /V <~ Ai:r d n = Ei so /(T d m p c 2 ). These two equations lead to 

N± 2T dlm m e c 2 fe^Y l/2 Y sc 



N 



rSC 



The flux ratio Y± c satisfies the following relation 



Y sc = y\c x _ Gj 



4 c (l+z) V^; 



r d f°° j dn ± 2 
fT / d 1± — 7 ± 



l+Y 



(A28) 



(A29) 



N± 1/2 1/2 P~ 1 ^ 

--> ml± 'Pp-2 ~ 
sc7±, c _ 2(T SC ) 2 

7c i+y i+y± ' 



SC 



:2y 



p-i y 
'p^2 i+y 



(A30) 



where y± = yj^+y]. , and we have used "f±x/lc = (l+Y) /(l +F±). The number ratio /V± /iV is equal to the ratio of the synchrotron 
fluxes of the pairs and the original electrons, F e± c /F Ec . Thus, the ratio of the eF £ fluxes is 

7 2 ±. p F £± A£±. m /£±,r l/2 (£±,p/£±, l nT 3/4 



llFeXeJec)- 1 ' 1 



ySC 

= l+Y± 

Equations (A27), (A30), and (A31) reduce to 

x(l — x— x 2 ) 



= -, Y sc = 



y| c = 2xy su , 



SC 



(A31) 



(A32) 



(l-2x-2x 2 ) 2 eg' ' l-2x-2x 2 ' 

which lead to a solution x ~ 0.34, y sc ~ 3.5, Y IC ~ 1.2, 7 ~ 4.7, y| c ~ 2.4, yf ~ 7.1, and Y± ~ 9.5. These factors are only 
dependent on and e e , being constant as long as 7± ;C < 7± >m < 7± )P < 7±,* is satisfied. 

Figure 3 shows an approximate spectrum of the emission by the original electrons and the pairs. The synchrotron spectrum of 
the pairs is approximated as 

te) te) fe) ■ 



/ \ V 4 / \ !/2 

(ifcj ' fOT £a < 6 < £± -- 



, for e ±;m < e < £ ±jJ ,, 
(^7) 4 " , for £±, p < £ < £±,t, 

l(3i)" W (*)" W . 



(A33) 



The SSC spectrum of the pairs is 



( ( > c \ 1/4 / \ 1/2 / \3 

{%) (*) fe)' for£<£ « 
1/4 , ,„ 

(*) 

1/4 



£F £ ^£ S ± c p F £ s ± C p X< 



for £ a < £ < £^ c m , 

(if;) 1 '". for 4 C m<£<4 c pI 
(jk-) 4 2 , for 4 c p < e < e 



(A34) 



e ±,p 
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£0.3 



-±Jfc! 



for £ > £^, 
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where e±~ m ~ 2j± . m e±, m and e*± k ~ 27^. t £±,i. The IC spectrum of the pairs is 

/ £ ic X 1 / 4 / xl/2, s3 
/ ic \ 1/4 / \l/2 

(^J j , for e a < e < s l £ m , 

forel m <e<el p , (A35) 
(if-;) 4+2 , for e% p < e < e 1 ^, 

/ ic \-f + 5 , N -£+l 

(&) (*) •*">■&■ 

where £± m ~ 2^± m e m and £± fc — 27j_ ££„,. The IC spectrum of the original electrons is the same of that of the pairs except for the 
flux normalization. We have superposed these two IC components in Figure 3. The energies e a and e 11 should be re-calculated by 
taking into account the emission of the pairs. For the adopted parameters, however, we found that e a is larger, and e~ n is smaller, 
than the original values only by a factor of ~ 2, which do not significantly affect the overall spectrum. We have e 77 > £± p , so 
that no further significant cascade emission is expected for the parameters adopted here. 

We have assumed that the energies of the pairs created in the emitting region are just given by the annihilated photon energies 
and the shock acceleration of the energetic pairs is not effective. The total pair energy density is comparable to the total injected 
energy density of the original electrons, while the ratio of the number densities is given by Equation (A29). Thus the characteristic 
Lorentz factor achieved by the pairs from their shock acceleration is estimated as r )±. ac — "/ m N/N± ~ 7±,/>(£ 77 / '£„ c ) 1 ^ 2 [(l + 
Y)/Y sc ] ~ J±, P - This implies that the shock acceleration does not affect significantly the pair energy distribution given from the 
annihilating photon energy distribution. 



-IC 



ic 
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